Effect of acid proticity on the thermodynamic parameters of charge transfer resistance in corrosion and passivation of nickel based glass alloy

The effect of temperature on electrochemical properties of Ni82.3Cr7Fe3Si4.5B3.2 glassy alloy in different acid proticity has been investigated utilizing AC and DC methods. Firstly, the handling of experimental data on the temperature dependence of charge transfer resistance, as well as corrosion current density permits us to determine the values of classical Arrhenius parameters as well as the thermodynamic ones considered approximately independent of temperature. This leads us to deduce a global interpretation on the phenomenon of corrosion and polarization. Secondly, the deviation to the linearity of the Arrhenius behavior and the real dependence on temperature of the thermodynamic parameters, permit us to clearly quantify the effect of the acid proticity and define, for the first time, the concept of current Arrhenius parameters and the current thermodynamic ones, as well as the modeling of the enthalpy–enthalpy compensation. Moreover, the effect of temperature can be investigated using the Vogel–Fulcher–Tammann model to reveal that the corresponding Vogel temperature has an interesting physical meaning.


Materials and methods
Ni-based alloy ingots Ni 82.3 Cr 7 Fe 3 Si 4.5 B 3.2 (wt%) alloy was prepared by rapid solidification supplied as ribbons of about 2.5-7.5 mm length and 20-50 µm thickness by Vacuumschmelze after polishing operation.
The electrochemical cell included three electrodes such as: the corroded sample as working electrode, a platinum-wire counter electrode and, a saturated calomel reference electrode for which each experiment is realized using a new alloy strip after degreasing in alcohol, rinsing with be-distilled water, and ultrasonic cleaning.
Because of the high aggressivity of acids in the studied temperature range, electrochemical characterization was performed in acidic solutions of 3.0 mol/L of HCl, H 2 SO 4 , or H 3 PO 4 for examining the electrochemical behavior of Ni-based glassy alloy at different temperatures ranging from 20 to 80° before the specimens were destroyed notably at higher temperature.
The electrochemical measurements, on samples with 2 cm 2 surface area, were performed by direct and alternating currents at several temperatures with an ACM Gill AC instrument.Measurements by electrochemical impedance spectroscopy were carried out at frequencies from 0.1 Hz to 30 kHz utilizing sinusoidal wave of 5 mV amplitude.Polarization measurements were realized, at a scan rate of 2 mV/s, from (− 800 to 2000) mV.
Complementary details are presented in our previous works 1, 22 .

Deviation to Arrhenius behavior Correlation between Arrhenius parameters
Activation energy control occurs when the electrode kinetics or corrosion rate is controlled by a slow electrochemical step.The activation energy (Ea) can be determined from Arrhenius plot that often used to analyze the effect of temperature on the rates of chemical reactions where the reactants molecules colliding probability becomes higher and the reaction proceeds faster at higher temperatures.In our study Ea of the corrosion process of glassy alloy in 3.0 mol/L solutions of the studied acids was obtained from the assumed linear variation of the corrosion reaction rate assigned as the reciprocal of the charge transfer resistance, 1/R ct , with temperature 1, 22 .The variation of the charge transfer resistance (R ct ) as a function of temperature (T) is considered following the Arrhenius-type equation (Eqs. 1 or 2).In this context, in previous works 1, 22 we have applied the Arrhenius-type equation in this exponential and logarithm form where the two Arrhenius parameters are generally supposed both constants practically independent of temperature.
where Ea is the activation energy, A ct is the pre-exponential factor (Table 1).Though the theories of collision, transition state, statistical physics, theory or chemical reaction rate have detailed and expressed the effect of temperature, generally 1, 22 , experimenters assumed that these two Arrhenius parameters are both constants and practically independent of temperature, where the plot of (lnR ct ) as a function of the reciprocal of absolute temperature (1/T) gives approximately a straight line 1,22 whose the slope is equal to (Ea/R) and the intercept on the ordinate is equal to (lnA ct ).Additionally, we can append that in addition to the y-intercept (lnA ct ) in (Fig. 1 of Ref. 22 ), we can also speak of the x-intercept (T A = -Ea/(R•lnAct)) previously named as the Arrhenius temperature 22 .In addition, we have shown that (− R•lnA ct ) is considered as an entropic factor (Fig. 3 of Ref. 22 ) and closely correlated with the activation entropy ΔS determined from impedance measurements in our earlier work 1 (Table S1).
(1) R ct = A ct e Ea RT (2) By similarity of some recent studies on the viscosity of liquid state which find that the Arrhenius temperature (T A ) is strongly correlated to the boiling temperature (T b ) [29][30][31] , and regarding that the BMGs are in solid state, we can assume that the Arrhenius temperature (Table 1) is probably in causal correlation with the corresponding glass transition temperature (T g ) which is strongly correlated with the melting point (T m ) 14,26,27 .This allows us in the future to make certain predictions, estimations and comparisons of different BMGs features.
Furthermore, the mutual correlation between the two Arrhenius parameters (Ea and lnA ct ) (Table 1), shows quasi-linearity inter-dependence expressed as follows: where (E a0 = 21.64 kJ mol −1 ) is the activation energy corresponding the zero value of the entropic factor, and the slope (τ 0 = 304.45K) is a temperature characteristic of the studied system at such conditions.Numerical applications require the use of the following convenient simple formulas: We point out that the values of coefficients can vary according to experimental conditions and to the system characteristics and specificity.
We notice that in the general case when the two Arrhenius parameters (Ea and lnA ct ) depend slightly on the temperature, the characteristic temperature (τ 0 ) becomes the derivative of the activation energy with respect of the entropic factor (Eq. 4) at constant pressure.
Since the observed linear dependence between the two Arrhenius parameters Ea and lnA ct (Eq.3), the charge transfer resistance R ct can be expressed only against of Ea or lnA ct (Eqs.5 and 6), respectively.Indeed, combining (Eqs. 2 and 3) we can write the following expressions:  The values of these coefficients can vary according to experimental conditions and to the characteristics and specificity of the studied system.In addition these two expressions are interesting when only one Arrhenius parameter is predicted by certain theory, the charge transfer resistance R ct can be then estimated by (Eqs.5 or 6).
To our knowledge, there is no theoretical and physical basis of this observed causal correlation or any developed estimative techniques for our original suppositions.Then, we will be able to give more closely our checking after applying these suggested equations by several researchers in the future.
Based on the studied correlation 26,27 between the glass transition temperature (T g ) and the melting temperature (T m(AxBy) ) of the base components of various metallic glasses where (300 K < T g < 900 K) and (600 K < T m(AxBy) < 1600 K) and the suggested expressions for the viscosity-temperature dependence [29][30][31] , we propose similar predictive expression form for the charge transfer resistance-temperature dependence expressed as follows.
As quick application of this formula, if we take (T m = 1726 K) for some Ni-MGs in (Fig. 1 of Ref. 26 ) and the Arrhenius temperature (T A = 517.59K) in the monoacid medium for our Ni-MG (Table 1), we find, via the (Eq.7), approximately a value of (T g = 739 K) which is included in the range of the studied set of Ni-MGs in (Fig. 1 of Ref. 26 ).We conclude that this approximative estimation can open an interesting path for future investigators to test, valid and improve the proposed model (Eq.7).
In the same way, by analogy with the viscosity-temperature dependence at liquid phase [29][30][31] which authors have discovered that the A ct is the pre-exponential factor (Eq. 1) is very close to the viscosity of the same system at vapor phase and at normal boiling temperature, we can presume that our pre-exponential factor (A ct ) in (Table 1) probably represents approximately the charge transfer resistance (R ct ) in liquid phase at the melting temperature (T m ) or at the glass transition temperature (T g ) under atmospheric pressure.

Effect of protons' number (x H )
For this query, and to test the effect of number of protons (x H ) of the acids (H x B), we plotted the two Arrhenius parameters against the logarithm of the acid proticity (lnx H ). Figure 1 shows a spectacular linearity, which leads us to wonder, for future investigations, if this relationship is valid for other strong acids with the same proticity.The slope of the straight line in Fig. 1a represents the activation energy gap ε g (Table 2) and it corresponds to the jump of energy value when the proticity increases with unity.The same ascertainment is valid for the Fig. 1b concerning the entropic factor gap σ g (Table 2).This observed linearity can be expressed as follows: where the straight line parameters are the intercepts on the ordinate (E a1 ) and (-R•lnA ct1 ) correspond to the values of the activation energy and entropic factor related to the monoacid HCl (x H = 1), respectively (Table 2).Numerical applications require the use of the following convenient simple formulas: 45  ) The values of these coefficients can vary according to the experimental conditions as well as the system characteristics and specificity.
However, Eqs. 1, 8 and 9 can be re-expressed to explicit the effect of protons number (x H ) of the polyacids (H x B) based on monoacid HCl (x H = 1) data.
where γ g represents a kind of free energy gap of acid protonation, R ct1 (T) is the charge transfer resistance related to the monoacid HCl (x H = 1) at given temperature (T) and (A ct1 ) is the corresponding pre-exponential factor (Eq. 1).
Numerical applications require the use of the following convenient simple formulae: The values of these coefficients can vary according to the experimental conditions as well as the system characteristics and specificity (Table S2).
Moreover, using linearization technique, we detect another interesting strong correlation between the Arrhenius temperature (T A ) and the compensation temperature (T comp ) with the acid proticity (x H ) depicted in Fig. 2.
Linearization technique leads to simple numerical applications requiring the use of the following convenient simple formulas: and (10)  where T comp =|ΔH°/ΔS°| And the six numerical values are all equivalent to absolute temperatures expressed in Kelvin.We hope that in future works we will find some further interpretations for probable physical significance of the amount of these values, especially for (T 0 ) indicated in the first member of the two previous equations (Table S3).

Activation energy-temperature dependence
Nevertheless, for the variation of the logarithm of charge transfer resistance (lnR ct ) as a function of the inverse of the reciprocal temperature (1/T), we have observed feeble net deviations from the linearity of the Arrhenius behavior in some results in literature 1, 22,33 .Generally, the linear regression used by the majority of experimenters is just an approximation, because the experimental points are not aligned in practically all researches.For this reason, we consider the slight variation of activation energy with the temperature around the constant values calculated by the classical Arrhenius-type equation, and then the interpretation of the sense of variation against temperature will be very interesting and enrich classical conclusions, discussions and explanations.To differentiate between the classical Arrhenius parameters independent of temperature and those of the new concept, we will name them as current Arrhenius parameters Ea(T) and lnA ct (T) which are dependent on temperature.So, Eq. 1 is dropped and replaced by a similar expression whose parameters become temperature dependent (Eq.13).
For that, we propose, as optimization by nonlinear regression, to simply fit the experimental results (lnR ct ) with respect of (1/T) with only a small-degree polynomial, which can be expressed in its general form as follows: In fact, in our situation, we are satisfied with two-degree polynomial (a 3 = 0, etc.) where results are given in Table 3 22 .
We note that, the mathematical derivation of Eq. 14 can lead us to determine (Ea) and (lnR ct ) using the following equations, In our case, the deviation to the linearity is feeble, so we will consider only a second degree polynomial in our nonlinear regression and consider with excellent approximation that the third coefficient is zero (a 3 = 0) 1 .Numerical applications require the use of the following convenient simple formulas: where we can inject, for each acid, the values of (a i ) from Table 3.The values of these coefficients can vary according to the experimental conditions as well as the system characteristics and specificity.
Results of nonlinear regression are given in Table 4. Figure 3 illustrates this interesting variation of Arrhenius parameters with temperature.We add that we will be forced to augment the polynomial-degree when the general trend of data points, in the plot of lnR ct = f(1/T) has a strong curvature or a change of curvature (inflection point).
Graphical analysis of (Fig. 3) shows that values of the current Arrhenius parameters Ea(T) and lnA ct (T) meet those of classical Arrhenius parameters Ea and lnA ct (Table 4) for practically a unique temperature (T cr = 320.71K,47.56 °C) named as crossover temperature which is approximately equal to the arithmetic mean (319.15K,46 °C) of the set of five working temperatures indicated in (Table 4).We can notice that this temperature can be an excellent working temperature for our studied system because it gives precise values of Arrhenius parameters whether we treat the data with linear or non-linear regression.
We add that the crossover temperature (T cr ) is not specific for the system like the characteristic temperature (τ 0 ) studied in previous work 1 , one must be careful because it cannot be characteristic of the studied system; it is only an intermediate mathematical variable obeying Eq. 19 by simply indicating approximately the temperature of the system during its process.The Eq. 19 can be easily obtained by equalizing Eq. 1 for classical parameters and Eq. 17 for current ones.
where (Ea) represents the activation energy obtained by linear regression.For numerical applications, we can inject, for each acid, the values of (a i ) from Table 3 and those of (Ea) obtained by linear regression from Table 4.Moreover, we can also use the values of energy parameters determined in our previous work 1 as follows, We notice that the middle-temperature (T md = 323.15K,50 °C) of the studied temperature range is very close to the temperature parameter (T cr ) as it is shown in Table 4. Indeed, to give an approximate estimation of the mean ( 19) Results of values of the two current Arrhenius parameters Ea(T) and lnA ct (T) obtained by Eqs. 17 and 18. (a): Crossover temperature T cr (Fig. 3).(b); current Arrhenius parameters Ea(T) and lnA ct (T) calculated with Eqs. 17 and 18, respectively.(c): logarithm of charge transfer resistance (lnR ct ) estimated with Eq. 13. (d): from 22 .(e): Standard deviation σ.Which can be adapted for Eq.21 and can lead to a convenient expression (Eq.22) for an average value of activation energy without using direct linear regression of lnR ct with 1/T.By similarity to the observed mutual correlation between the two Arrhenius parameters (Ea and lnA ct ) 22 , likewise we have thought about inspecting the mutual dependence between the two current Arrhenius Ea(T) and lnA ct (T) by plotting one parameter against the second for the three studied polyacids H x B. In fact, the Fig. 4 shows interesting causal correlation whether for each acid separately or for all three together, for which the quasi-linearity inter-dependence can be expressed as follows: where E a0 (Table 5) is a current activation energy corresponding the zero value of the entropic factor, and the slope τ 0 (Fig. 4) is equivalent to an absolute temperature characteristic of the studied system under such conditions and within the working temperature range.Hence, some experimenters 1 interpret the sign of the deviation to the Arrhenius linearity as a sub-Arrhenius or super-Arrhenius behaviors.Therefore, we notice that the activation energy Ea (Eqs.15 and 17) can be interpreted as a potential energy barrier which is assumed to be dependent on temperature Ea(T).
In case of positive values of the activation energy derivative with respect to the reciprocal of absolute temperature at constant pressure: ∂E a (T) ∂(1/T) P , super-Arrhenius behavior is observed, whereas for negative values, it is the sub-Arrhenius behavior.For this derivation, the values neighboring zero lead to the classical Arrhenius behavior (Eq.17) and the a 1 -value tends to the classical Arrhenius activation energy Ea, independent of temperature (Table 4).
Or else we can write the following: In our case, the system exhibits a super-Arrhenius behavior, i.e. the potential energy barrier is reduced whenever the temperature rises.

Vogel-Fulcher-Tammann model
Generally, in case of clear deviation to the Arrhenius behavior, experimenters try to classify their results in the super-Arrhenius behavior or the sub-Arrhenius one 22 or apply the Vogel-Fulcher-Tammann-type equation (VTF or VFT) 14,[23][24][25][26][27][28] especially when there is some divergence of experimental values at low temperature.
However, the variation of the logarithm of charge transfer resistance (lnR ct ) against the inverse of the reciprocal temperature (1/T) exhibits a feeble deviation from the linearity of the Arrhenius behavior in our previous work (Fig. 1 of 22 ) and in some results in literature 26,[32][33][34] .In addition, regarding the suggested pseudo-hyperbolic behavior in previous work 22 where the divergence of the variation of charge transfer resistance (R ct ) for low temperatures (Fig. 6 of 22 ) exhibits a kind of vertical asymptote, we propose to explore the Vogel-Fulcher-Tammanntype equation (VTF) which is characterized by a vertical asymptote and is used when the test of Arrhenius behavior fails 14,[23][24][25][26][27][28] .In the case of the nonlinear behavior, it is found that the temperature dependence of charge transfer resistance can be physically fitted with the frequently VTF-type equation 14,[23][24][25][26][27][28] expressed as follows: where A 0 and B 0 are optimal constants and T 0 is the Vogel temperature.It's also interesting to use the modified VTF equation which is expressed as follow: where R is the perfect gas constant, E 0 is the VTF activation energy and, A 0 and, T 0 , are the pre-exponential factor and the Vogel temperature generally comparable to the glass transition temperature in viscosity property 14,[23][24][25][26][27][28] .Experimental data are presented in Table 6 and presented in Fig. 5.For numerical applications, we can inject in Eq. 27, for each acid, the values of (T 0 ), (lnA 0 ) and (E 0 ) from Table 6.
The Eq. 27 shows that when the temperature (T) tends toward (T 0 ), it implies that the charge transfer resistance (R ct ) becomes infinity, which indicates that the corrosion is inhibited.We see that (T 0 ) decreases whenever the proticity (x H ) of the acid increases (Table 6) while the VTF-energy (E 0 ), which is in close relation with the activation energy (Ea), varies in the reverse sense.

Correspondence between the two models
Generally, several experimenters manipulate physical and chemical quantities using various models for comparison and to develop discussions, interpretations and conclusions.This section falls within the scope of the correspondence between the modified Arrhenius equation (Eq. 3) and the VTF model (Eq.27) to optimize the number of the models commonly utilized for investigations.
When we inject the second member of the VTF expression (Eq.27) in Eqs. 15 and 16 which are deduced from the principal modified Arrhenius equation (Eq.2) we can find the new expressions (Eqs.14 and 15) of the Arrhenius parameters (Ea) and (lnA ct ) by comparing term by term without the direct use of the classical Arrhenius-type equation (Eq.2).
Analyzing the expression of (Eqs.14 and 15) we can conclude that the two parameters (E 0 ) and (lnA 0 ) of the VTF model are simply mathematical intermediary of calculation, except the third one (T 0 ) which has physical significance indicating that the charge transfer resistance value (R ct ) diverges and reaches a very high value when the temperature (T) is very close to (T 0 ).Moreover, we can say that the non-zero (T 0 )-value is the principal cause of the deviation to the linearity of Arrhenius behavior.In fact, we can see in (Eqs.14 and 15) that when (T 0 ) tends to zero, the two parameters (E 0 ) and (lnA 0 ) become identical to those of Arrhenius (Ea and lnA ct ) and in only this situation the two VTF parameters (E 0 ) and (lnA 0 ) have a full physical meaning.Furthermore, some theorists state that the Vogel temperature (T 0 ) is in causal correlation with the corresponding glass transition temperature (T g ) 14,[23][24][25][26][27][28] .

Thermodynamic parameters-temperature dependence
Treatments of the free Gibbs energy-temperature dependence considering the corresponding thermodynamic parameters (ΔH°) and (ΔS°) are rarely done in literature [32][33][34] .The majority consider approximately the constancy of theses parameters and interpret globally the phenomenon governing their studied systems based on their signs (positive or negative), their amounts (high or low values) and not on the eventual slow variation with temperature.
Then, regarding the Eq.32, the Eq. 31 can be reformulated as follows: So, the plot of R•ln(N A hR ct /RT) as a function of 1/T gives a quasi-straight line, with directly a slope of (ΔH°) and an intercept to the ordinate of (-ΔS°).Results are presented in Table 7 and depicted in Fig. 6.The negative values of the Gibbs free energy (ΔG°) indicate that the process of the intermediate complex in the transition

Effect of protons' number (x H )
Similarly with our previous work 1, 22 for the Arrhenius parameters, we will test the effect of the protons number (x H ) of the polyacids (H x B).For that, we plotted the two thermodynamic parameters (ΔH°) and (ΔS°) against the logarithm of the acid proticity (lnx H ). Figure 7 shows a linearity, which leads us to wonder if this relationship is valid for other strong acids with same proticity.The slope of the straight line in Fig. 7 represents the enthalpy gap ε g (Table 8) and it corresponds to the jump of energy value when the proticity increases with unity.The same ascertainment is valid for the Fig. 7 concerning the entropy gap σ g (Table 8).This observed linearity is similar to the Eqs.8 and 9 with the same parameters values (Table 2) and it's expressed as follows: where the straight line parameters (ΔH°1) and (ΔS°1) correspond to the values of the enthalpy activation and entropy activation related to the monoacid HCl (x H = 1), respectively (Table 8).Numerical applications require the use of the following convenient simple formulas: We expect that the values of these coefficients can vary according to the experimental conditions as well as the system characteristics and specificity.
However, Eqs. 30, 34 and 35 can be re-expressed to explicit the effect of number of protons (x H ) of the polyacids (H x B) based on monoacid HCl (x H = 1) data.
where γ g represents the free energy gap of acid protonation (i.e. when proticity has changed by one unit), ΔH°1 is the enthalpy of activation related to the monoacid HCl (x H = 1) at given temperature and ΔS°1 is the entropy of activation.Numerical applications require the use of the following convenient simple formulas: The values of these coefficients can vary according to experimental conditions as well as the system characteristics and specificity (Tables S1, S2 and S3).
We expect that these suggested expressions can be utilized in future estimations for other experimental conditions or other studied materials.We conclude that the proticity (x H ) of the polyacid (H x B) has a substantial effect, which can be modeled for future prediction, or estimation and can induce theorists to develop or improve theories.

Correlation between the Arrhenius parameters and thermodynamic ones
Analysis of the enthalpy of activation (ΔH°)-values and those of the Ea, in the (Fig. 8), shows that the Ea and ΔH° values are very closely related.The same conclusion is also attributed to the correlation between the Arrhenius entropic factor (-R•lnAs) and the entropy of activation (ΔS°).Starting from the fact that the values of the two slopes of (Fig. 8) are practically equal to unity (1.00004 and 0.999994), the following expressions are proposed.
where (δH° = 2.675 kJ mol −1 ) and (δS° = 330.46J K −1 mol −1 ) are the enthalpy increment and the entropy increment respectively.Numerical applications require the use of the following convenient simple formulae: You have to be careful during discussions, interpretations and comparisons between amount values obtained for (R ct /Ω cm 2 ) or (R ct /Ω m 2 ), because the use of CGS or SI systems, which differ in the scale of base units, during the calculations of Arrhenius parameters or thermodynamic ones, doesn't affect the values of the activation energy (Ea) and the activation enthalpy (ΔH°), while it gives difference of (± R•ln10 4 = ± 76.579 J K −1 mol −1 ) when calculating the Arrhenius entropic factor (-R•lnA ct ) and the entropy of activation (ΔS°), and this is due to the conversion (cm 2 ↔ m 2 ).
We conclude that we can estimate one parameter when the other one is determined by any other technique.We notice that this shift is also due to the fact that the Arrhenius parameters represent the movement between two energy levels related to transition states, while the thermodynamic parameters, as state functions, represent www.nature.com/scientificreports/ the movement between two energy levels related to equilibrium thermodynamic states 1,22,30-37 .We add that the slopes values of the two straight lines of Fig. 8 are very near to the unity explaining then why that (Ea and ΔH°) and (-R•lnA ct and ΔS°) have approximately the same value of gap or "jump" (ε g ) and (σ g ) when the number of protons (x) of the acid changes by one unity (Tables 2 and 8).Furthermore, mathematical comparison considering the expressions of Eqs. 1 and 2 and the Eqs.30-33, which partly include both terms of Eqs. 1 and 2 in each member leads us to expect that the enthalpy increment (δH°) is in close correlation with the contribution of thermal agitation on the activation enthalpy of the thermal stability related to the spontaneous formation of the hard non-reactive surface of passive film that inhibits the further corrosion (Tables S1, S2 and S3).

Mutual correlation between the thermodynamic parameters
Analysis of the mutual correlation between the enthalpy of activation (ΔH°) and the entropy of activation (ΔS°) represented by the (Fig. 9) shows an excellent linearity which can permit us to estimate the free energy (Eq.32) with only one thermodynamic parameter (ΔH°) or (ΔS°).This issue is interesting when we can theoretically predict one parameter; we can then deduce the value of the other one by this observed linearity (Tables S1, S2 and S3).So, the linear enthalpy-entropy dependence observed in Fig. 9 can be expressed as follows, where (τ 0 = 304.45K) and (∆S o c = 250.618J K −1 mol −1 ) are the characteristic temperature (Eq. 3) and the characteristic entropy of the studied system at such conditions, respectively.We can consider that the (τ 0 ∆S o c )-product (Eq.42) is equivalent to a characteristic enthalpy (∆H o c = 76.301kJ mol −1 ), the Eq.41 can be rewritten as follows, Numerical applications require the use of the following convenient simple formulas: We deduce then, each investigated system has two main specific independent parameters (∆H o c ) and (∆S o c ), and a dependent parameter (τ 0 ) which can be simply deduced by the Eq.42.We note that (∆S°c) corresponds theoretically to the limit of the endothermicity, i.e. the activation entropy for which the activation enthalpy becomes zero and changes sign.

Case of temperature-dependent thermodynamic parameters
Almost, all researchers fit thermodynamic behaviors of experimental data in linear regression to conclude about the global thermal character of the studied process (i.e.endothermic, etc.).Analyzing the abovementioned study of Arrhenius parameters-temperature dependence, we conclude that is better if we think about the fitting by nonlinear regression to reduce the discrepancy between the straight line and the uncertainty bars of some of experimental scatter points in order to obtain more accurate values of the two thermodynamic parameters (ΔH°) and (ΔS°).Then, the linearity of (ΔG°) with the absolute temperature (T) of Eq. 33 and that of (ΔG°/T) with the reciprocal of absolute temperature (1/T) are abandoned and replaced by a polynomial equation with two or three degrees expressed as follows, Given the deviation to the linearity of (ΔG°/T) with the reciprocal of absolute temperature (1/T) observed in literature [29][30][31] is generally not very significant, we judge that we need to consider only a two-degree polynomial fitted in our nonlinear regression (Fig. 10) and acquire excellent approximation where the third coefficient is zero (a 3 = 0) 1, 22 .
The corresponding a i -coefficients of nonlinear fit are presented in the Table 9.We can clearly see that the statistical quality has improved quite well, and despite this we must pay attention to the fact that this could also be due to the small number of used working temperatures.
The two main thermodynamic parameters, such as the enthalpy (ΔH°) and the entropy (ΔS°) can be determined from the basic thermodynamic Gibbs free energy relationship (Eq.32).
Te determine values of the two thermodynamic parameters (ΔH°) and (ΔS°), we have to plot the ratio Gibbs free energy by temperature (ΔG°/T) with the reciprocal of absolute temperature (1/T).The mathematical handling of Eq. 32 can lead us to determine values of the two thermodynamic constant parameters (ΔH°) and (ΔS°) using ( 42) the Eq.45 in the case of linear behavior (a 2 = 0 and a 3 = 0 in Eq. 44) where the slope value of the straight line is (ΔH°) and intercept on the ordinate is (-ΔS°).
Generally, the two thermodynamic parameters become variable with temperature and mathematical derivation of Eq. 32 can lead us to determine values of the two thermodynamic non-constant parameters ΔH°(T) and ΔS°(T) using the following equations.Application of Eqs.46 and 47 leads to the general following convenient general expressions.
Numerical applications require the use of the following convenient simple formulas: where we can inject, for each acid, the values of (a i ) from Table 9.The values of these coefficients can vary according to the experimental conditions as well as the system characteristics and specificity.Values of the two thermodynamic parameters (ΔH°) and (ΔS°) obtained by linear and nonlinear regressions are presented by Table 7 and depicted in Fig. 11 for the three acids.
Graphical analysis of (Fig. 11) shows that values of the two current thermodynamic parameters ΔH°(T) and ΔS°(T) meet those of classical thermodynamic parameters (ΔH°) and (ΔS°) (Table 7) for practically the same unique temperature (T cr = 320.70K,47.55 °C) named as crossover temperature which is approximately equal to the arithmetic mean (319.15K, 46 °C) of the set of five working temperatures indicated in (Table 7).We can notice that this temperature can be an excellent working temperature for our studied system because it gives precise values of thermodynamic parameters whether we treat the data with linear or with non-linear regression.
We add that the crossover temperature (T cr ) is not specific for the system like the characteristic temperature (τ 0 ) studied in previous work 1 , one must be careful because it cannot be characteristic of the studied system; it is only an intermediate mathematical variable obeying Eq.50 by simply indicating approximately the temperature of the system during its processing.The Eq. 50 can be obtained easily by equalizing Eq. 45 for classical parameters and Eq.48 for current ones.
where (ΔH°) represents the enthalpy of activation obtained by linear regression (Table 7).
For numerical applications, we can inject, for each acid, the values of (a i ) from Table 9 and those of (ΔH°) obtained by linear regression from Table 7.We notice that the middle-temperature (T md = 323.15K;50°C) of the studied temperature range is very close to the temperature parameter (T cr ) as it is shown in Table 7. Indeed, to give an approximate estimation of the mean activation enthalpy ΔH° that should be obtained by classical linear regression we can apply the following reasoning.We can calculate the average value of the function ΔH°(T) expressed by Eq.48 over a temperature domain [T min ,T max ] such as in our situation; [293.15,353.15]K,by the following expression: Which can be adapted, using Eq.48, and can lead to a convenient expression (Eq.52) for an average value of enthalpy of activation without using direct linear regression of (ΔG°/T) with 1/T.Inspired by to the observed mutual correlation between classical thermodynamic parameters (ΔH°) and (ΔS°) abovementioned, we similarly thought of inspecting their mutual dependence between the two current thermodynamic parameters.ΔH°(T) and ΔS°(T) by plotting one parameter against the second for the three studied polyacids H x B. In fact, the Fig. 12 shows interesting causal correlation, analogous to that of Eq. 43, (51) www.nature.com/scientificreports/whether for each acid separately or for all three together, for which the quasi-linearity inter-dependence can be expressed as follows: where ∆H o c (Table 10) is a current activation enthalpy corresponding the null value of the activation entropy, and the slope τ 0 (Fig. 12) is equivalent to an absolute temperature characteristic of the studied system at such conditions.
Another advantage, from which we benefit from the dependence on temperature, and other than that of the deduction of the endothermic character of the process, is that which indicates the sense of variation of the enthalpy with the temperature (Table 7) to know if this endothermicity is accentuated or attenuated and thus allowing us to choose the optimal working temperature.The same goes for the discussion and interpretation of entropy.In other words, the dependence with the temperature offers more explanations and fine interpretations than those deduced from only the constancy of thermodynamic parameters.

Conclusion
In one previous work 1 , we have studied the effect of the temperature as well as the nature of polyacid medium (HCl, H 2 SO 4 and H 3 PO 4 ) on the tendency toward passivation and corrosion resistance for the nickel based glassy Ni 82.3 Cr 7 Fe 3 Si 4.5 B 3.2 alloy.
In another previous work 22 , we have taken into account the deviation to the linearity of the Arrhenius behavior and modeled this particular dependence on temperature.Consequently, we have ended up at the following points.(i): The mutual correlation between the Arrhenius parameters, such as, the activation energy determined from polarization and impedance measurements and the entropic factor of Arrhenius permit to reveal a new concept of the current Arrhenius temperature.(ii): The causal correlation between the Arrhenius parameters and thermodynamic parameters, permit us to conclude that the activation energy can be considered, with a reliable approximation, as a thermodynamic function.(iii): The deviation to the Arrhenius linearity can be classified as a super-Arrhenius behavior for the studied system.(iv): The effect of the proton number in polyacid is quantified and modeled.(v): The novel suggested model describing a pseudo-hyperbolic behavior gives an excellent agreement (R-square ≈ 1) with the experimental data for each acid separately or for all three together.
As continuation of the abovementioned investigations, we keep on presenting furthermore original modeling of the temperature effect on the thermodynamic parameters as well as the effect of the polyacid proticity.Then, in the present work we extricate the following key points.(i): The mutual correlation between the Arrhenius parameters, such as, the activation energy determined from polarization and impedance measurements and the entropic factor of Arrhenius permit to rewrite the Arrhenius-type equation with only one parameter Ea or lnA ct (Eqs.5 and 6) and to facilitate, for theorists, to predict the value of one parameter (Eq. 3) when the other one can be predicted by certain theory or approximation.(ii): The effect of the proton number (x H ) in polyacid (H x B) on the two Arrhenius parameters (Ea and lnA ct ) is modeled by original expressions (Eqs.8 and 9) for which we introduce the new concept of the activation energy gap ε g and the entropic factor gap σ g permitting to estimate new values of Arrhenius parameter Ea or lnA ct of a polyacid, using a power law expression (Eqs.10 and 11), when the same parameter of another acid is available (Table S2).(iii): The deviation to the Arrhenius linearity is simply modeled by a polynomial expression (Eq.14) which reveals the new concept of current Arrhenius parameters depending on temperature (Eqs.17 and 18) and then permits us to deepen the interpretation and discussion, and add additional elucidations of the effect of temperature on the Arrhenius parameters and the related phenomena governing the studied system.(iv): In the same context, we introduce the notion of crossover temperature (Eq.19)  as an optimal working temperature.(v): The mutual correlation between the current Arrhenius parameters, depending on temperature, exhibits the same behavior and interdependence whether for each acid separately or for all three together (Eq.23).(vi): We introduced for the first time the Vogel-Fulcher-Tammann-type equation to model the variation of the charge transfer resistance with the respect of temperature (Eq.27) which shows an excellent agreement with experimental data.(vii): In the same context, we have given expressions (Eqs.(1/T) permits us to express the dependence of the current thermodynamic parameters ΔH°(T) and ΔS°(T) on temperature (Eqs.48 and 49) and to propose the notion of crossover temperature (Eq.50) as an optimal working temperature.We conclude that even for the case of the current thermodynamic parameters depending on temperature, the interdependence maintains the linear behavior.(xii): we proposed predictive expression form for the charge transfer resistance-temperature dependence to approximately estimate the glass transition temperature (T g ).After all, it should be mentioned that the previous interpretations made based on the order-disorder or enthalpy-entropy compensation effect could be enhanced by raising the number of working temperatures.Likewise, it cans explain the tendency of these alloys to undergo a disorder-to-order transformation in certain temperature range.We note that some interfaces can introduce an order-disorder transition of this two-dimensional layered network into molecules, leading to increased diffusional characteristics and reduced bonding.Moreover, we can observe an order and disorder combined corrosion morphology of dual-phase Ni-based alloy in the passive state 22,[32][33][34][35][36][37][45][46][47][48][49] .
Finally, we must be wary, an interpretation with few experimental data may lead to less certain or unconvincing conclusion.

Figure 1 .
Figure 1.Variation against the logarithm of the acid proticity (lnx H ) of (a): the activation energy (Ea) and (b): the entropic factor of Arrhenius-R•ln(A ct ).

Figure 2 .
Figure 2. Variation against the acid proticity (x H ) correlation between (a): the Arrhenius temperature (T A ) and (b): the temperature of compensation (T comp ).

Figure 5 .
Figure 5. Variation of the logarithm of charge transfer resistance (lnR ct ) against 1/(T-T 0 ) for different acids.

28 and 29 )Table 10 .
permitting to link the VTF parameters and those of the Arrhenius-type equation to facilitate for the users to obtain double results by using only one chosen model.(viii): The effect of the proton number (x H ) in polyacid (H x B) on the two thermodynamic parameters (ΔH° and ΔS°) is modeled by satisfying linear expressions (Eqs.33 and 34) which permit us to estimate the parameters values to one polyacid, knowing those of the monoacid or another polyacid.(ix): The correlation between the Arrhenius parameters and the corresponding thermodynamic (53) �H o (T) = τ 0 • �S o (T) − �H c o Optimal thermodynamic parameters (∆H o c , ∆S o c and τ 0 ), from linear and nonlinear regression of Eqs.41 and 43.*Fit in linear regression for thermodynamic parameters as independent of temperature.**Fit in nonlinear regression for thermodynamic parameters dependent on temperature.
ones, such as (Ea and ΔH°) and (-R•lnAs and ΔS°), leads us to determine the value of one parameter, simply by a gap (Eqs.39 and 40) named enthalpy or entropy increment.(x): The mutual correlation between the thermodynamic parameters (ΔH° and ΔS°) shows a linear dependence (Eqs.41-44) and reveals values of characteristic enthalpy and characteristic entropy, specific to the studied system in such conditions.(xi): The feeble deviation to the linearity of (ΔG°/T) vs.

Table 1 .
Arrhenius parameters (Ea and lnA ct ), the entropic factor of Arrhenius-R•ln(A ct ) from linear regression of Eq. 2, the activation entropy ΔS and Arrhenius temperature (T A ). *Calculated via the variation of the charge transfer resistance (R ct /Ω cm 2 ) as a function of temperature (T).

Table 7 .
Resultsstate for the corrosion of Ni 82.3 Cr 7 Fe 3 Si 4.5 B 3.2 alloy in acidic medium is spontaneous, while the positive values of enthalpy (ΔH°) show that the process is endothermic 1,22,30-37 .

Table 8 .
Straight lines parameters from linear regression of Eqs.34 and 35.*Correspond to HCl experimental data, see Table 6.

Table 9 .
Optimal a i -coefficients values obtained by non linear regression for Eq.44 with second polynomial degree for glassy Ni 82.3 Cr 7 Fe 3 Si 4.5 B 3.2 alloy corrosion for the three polyacids HCl, H 2 SO 4 and H 3 PO 4 .